Phase-separated ferromagnetism in spin-imbalanced Fermi atoms loaded on an optical 
ladder: a density matrix renormalization group study 
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We consider repulsively-interacting cold fermionic atoms loaded on an optical ladder lattice in 
a trapping potential. The density-matrix renormalization-group method is used to numerically 
calculate the ground state for systematically varied values of interaction U and spin imbalance 
p in the Hubbard model on the ladder. The system exhibits rich structures, where a fully spin- 
polarized phase, spatially separated from other domains in the trapping potential, appears for large 
enough U and p. The phase-separated ferromagnetism can be captured as a real-space image of the 
energy gap between the ferromagnetic and other states arising from a combined effect of Nagaoka's 
ferromagnetism extended to the ladder and the density dependence of the energy separation between 
competing states. We also predict how to maximize the ferromagnetic region. 

PACS numbers: 03.75.Ss, 67.85.Lm, 71.10.Fd 



Why cold atoms on a ladder? — Ultracold atom sys- 
tems offer ideal opportunities for systematic studies of 
novel quantum many-body phenomena, since they are 
not only extremely clean, but also controllable to an un- 
usual degree [l|. Particularly, the inter-atomic interac- 
tion tunable with the Feshbach resonance create a unique 
stage for exploring strongly-correlated systems. If we fur- 
ther turn to cold atoms on optical lattices (OL's) pre- 
pared by standing waves from laser beams, they provide 
an even more versatile playing ground 0, 0]. Among 
the attractive targets are the metal-insulator transition, 
d-wave superfluidity, and various magnetisms. Actually, 
the antiferromagnetism via the superexchange interac- 
tions between localized Bose atoms on an optical super- 
lattice [3| and the Mott insulator in fermionic atoms on 
cubic OL [3, [1| have already been achieved. 

Now, one big issue in the field of strongly correlated 
systems is the itinerant ferromagnetism, which is, de- 
spite the long history, still far from fully understood. 
Cold atom systems, with their tunability, should be an 
unprecedented place for realizing the itinerant ferromag- 
netism. While there are many theoretical studies on itin- 
erant magnetism in ultracold atom systems [6MTT| . Jo et 
al. recently reported that the Stoner instability was ob- 
served in an ultracold fermionic atom system without OL 
12]. However, it has been pointed out that some factors 
other than Stoner instability may also be involved [l3j . 

In this Rapid Communication, we study strongly- 
correlated fermionic atoms with a spin-imbalance loaded 
on an optical ladder, and propose a way to realize and 
observe the ferromagnetism originated by the finite hole- 
density Nagaoka mechanism. One tunability unique to 
cold atoms is we can control the spin balance, which en- 
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FIG. 1. (Color online) Two-leg ladder optical lattice (corru- 
gated surface in blue) is schematically shown along with the 
atomic wavefunction (pale cloud) and spins (arrows). Inset 
depicts simple lattice (along x) and superlattice (y) poten- 
tials for creating the optical ladder. 



riches the quantum states. For example, a phase separa- 
tion between superfluid and normal phases was observed 
in attractively interacting ultracold Fermi atoms [HI, EH ■ 
On the other hand, a magnetic structure with spin im- 
balance in repulsively interacting fermionic atoms with 
and without OLs has been studied in a weak-correlation 
regime 

The reason why we take the ladder is as follows. The 
long history of itinerant ferromagnetism has made us re- 
alize that Stoner's picture is mean-field theoretic, while 
almost the only model in which itinerant ferromagnetism 
can be rigorously shown is Nagaoka's ferromagnetism, 
which requires rather a pathological (infinitely strong 
interaction and an infinitesimal doping in a half-filled 
band) limit [Hj]. The Nagaoka ferromagnetism has yet 
to be experimentally confirmed, although von Stecher et 
al. recently proposes a way to realize the Nagaoka fer- 
romagnetism in optical plaquette systems with the high 
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FIG. 2. (Color online) The particle density (n^; blue solid line) and spin density (si\ red dashed line) for the spin polarization 
p = (a) , 0.4 (b) , 0.6 (c), and 0.8 (d) with U = 50, V = 3, i± = 1, and L = 100. P, PM, and FF denote paramagnetic, 
partially magnetic, and fully ferromagnetic states, respectively. The dash-dotted lines represent the m at numerically obtained 
phase boundaries. 



tunability of OL system [11] . On the other hand, some 
theories indicate that the two-leg ladder can realize an 
itinerant ferromagnetism that accommodates finite inter- 
action and finite doping, which is connected to Nagaoka's 
in the limit of U — > oo or large inter- leg transfer [T3, [l8[ ■ 
Experimentally, an optical ladder may be created by a 
superposition of normal lattice potential (along x) and 
super-lattice potential (y) (Fig. [T]). This should be real- 
izable because a superlattice potential has been used in 
the experiments with bosonic atoms 0, [H| . To confine 
the atoms we impose a trapping potential as well, and we 
explore the ground states in this situation. We want to 
treat very strong repulsive interactions, but the system, 
being quasi-one dimensional, can be treated with the 
density-matrix renormalization-group (DMRG) method 
[2fl |2ll . which can also treat effects of the trapping poten- 
tial [22| . We show that an itinerant ferromagnetism does 
appear in an optical ladder with an interesting phase sep- 
aration between a fully-polarized ferromagnetic domain 
originating from the extended Nagaoka mechanism with 
finite hole density. This ground state should be easily 
observed by in situ imaging method [14j . which is the 
main reason why we consider the spin-imbalance. In ad- 
dition, we find the enhancement of the formation of the 
finite hole-density Nagaoka ferromagnetism by the spin- 
imbalance. 

Formulation — We take the simplest possible repul- 
sive Hubbard model on a ladder in a trapping potential 
applied along the leg. The Hamiltonian then reads, in 
the standard notation, 

a = -t - *-l Yl 

+ [/^n? t n£ + ^(i-ag 2 n?, (1) 

i,a i,a,a 

where, i, j = 1,2, ■ ■ ■ , L (L: the ladder length) label the 
sites along the legs with (i,j) being nearest- neighbor sites 
connected by a hopping integral t, a — 1,2 labels the sites 



along the rungs with t± being the hopping along the rung. 
U(> 0) is the repulsive, on-site interaction (23|, while V is 
the strength of the trapping potential, which is assumed 
to be harmonic around the center, x c = (L+l)/2 with the 
lattice constant taken to be unity. It should be mentioned 
that the maximum value of U for stable Hubbard model 
in the current experiment is 180 [4]. We assume that 
the harmonic confinement along the rung (superlattice) 
direction is very weak. In addition, we also assume that 
the inter-ladder hoppings are negligible. 

The total spin imbalance is defined as p = (Nf — 
NJ/N, where N a = and N = 5Z CT N a , and 

(•) represents the expectation value. All the numerical 
results respect the symmetry, (n\ a ) = (n^), so that we 
can introduce the rung particle density n, L = (nL) + (h}i), 
and the rung spin density Sj = (n}+) — ( n li)- For con- 
venience, we use dimensionless parameters as i± = t±/t, 
U = U/t, and V = V/t x 10 3 . In this paper, we fix the 
total number of atoms as N = 100 with the length of 
the ladder 80 < L < 180. The number, m, of the states 
retained in the present DMRG calculations is 800-1200, 
which are numerically shown to give converged results. 

DMRG results — Let us show the DMRG results. Fig- 
ure [5] shows profiles of m and Sj in the ground states for 
various p in a strongly-interacting case (U = 50) with 
t ± = 1 and V = 3. In the spin-balanced case (p = 0) 
[Figj2ja)], a paramagnetic (P) phase appears. As we in- 
crease p to 0.4, a partially magnetic (PM) phase starts 
to appear in the middle of the system [Fig. HJb)]. If we 
further increase p to 0.6, fully ferromagnetic (FF) regions 
emerge in the ground state [Fig. EJc)] in a more complex 
phase-separated structure. For p = 0.8, the FF region 
covers almost the whole system, sandwiched by narrow 
PM and P regions [Fig.^d)]. 

In order to demonstrate both of the effect of strong in- 
teraction and the effect of the ladder configuration (a con- 
nectivity condition [16]) are at work for realization of the 
FF phases in the regions with 0.8 < n < 0.9, we compare 
the result with the one in a weakly-interacting case with 
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FIG. 3. (Color online) The particle density m (blue solid 
lines) and spin density s, (red dashed lines) for a weakly- 
interacting ladder system (U = 1, V = 2, tx = 1, and L = 
100) with p — 0.6 (a) or 0.8 (b), and for a strongly-interacting 
single chain (U = 50 and V = 2.5) with p = 0.6 (c) or 0.8 (d) 
for N = 50. 
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FIG. 4. (Color online) The density dependence of the energy 
gap AE = (_Es.=o — Es, = N/2)/Lt in homogeneous systems 
(V = 0) with tx = 1 for various values of (7. The inset is 
the same plot for the homogeneous single chain. Here the 
parameters are L — 50, tx = 1 and m = 1200 (ladder) or 
m = 1000 (chain). 



connectivity [two- leg ladder, U — 1, V — 2, and ij_ = 1; 
Figs. EIa)(b)], and with the one in a strong- interacting 
case without connectivity [single chain, U = 50, and 
V = 2.5; Figs. Etc) (d)]. We can immediately see that FF 
phases do not emerge in these systems, although there 
are narrow FF regions around the edges. This clearly 
indicates that the FF phases with 0.8 < n < 0.9 in the 
ladder arise due to a combined effect of the strong in- 
teraction and the connectivity condition, i.e., the finite 
hole-density Nagaoka ferromagnetism [17| ■ 

Now, the question is what determines the spatial dis- 
tribution of majority and minority spins as the spin im- 
balance is increased. As a key, we can examine the energy 
gap, AE = (E Sz= o — E s ^ =N / 2 )/Lt, of the FF above the 
lowest energy of the nonmagnetic state in the uniform 
system (V = 0), as in the t-J model in [l?} but here 
for the Hubbard model. The gap depends on the par- 
ticle density n as depicted in Fig. As U is increased 
the curve approaches the AE = axis around n ~ 0.8, 
which means that the ferromagnetic phase becomes closer 
to the ground state. Such a non-trivial structure is natu- 
rally absent in the single-chain Hubbard model, as shown 
in the inset of the figure. 

With this picture we can interpret the results shown 
in Figs. [21 (21[c) and[2Kd). When the trapping potential 
is not too strong so that the spatial variation of the den- 
sity is slow, we can locally define the energy gap as given 
in FigUl The energy gap and the trapping potential de- 
termine the phase separation: Let us first look at the 
strongly interacting ID chain, both of the majority (= 
{rii + Sj)/2) and minority (= (rii — Si)/2) spins are accu- 
mulated in the middle where n ~ 1 [Figs. (3][c)(d)]. Non- 
FF states are preferred because both of AE(n) (Fig. 21 
inset) and the trapping potential favor n ~ 1, which is a 
realization of the Lieb-Mattis theorem (24[. This leads 
us to a natural interpretation that the ground-state en- 
ergy is lowered by gathering minority spins around the 



trap center to construct an antiferromagnetic correlation. 
On the other hand, in a strongly-interacting ladder, we 
obtain the ferromagnetic region, which may seem coun- 
terintuitive, but is in fact viewed as a manifestation of the 
extended Nagaoka mechanism [l7| . Namely, the minor- 
ity spins in Fig. [U[d) reside in the regions with relatively 
high potential energy (i < 31 and 70 < i), but this is 
compensated by a large energy AE gained by forming 
the non-FF state (Fig. HJ) . Conversely, FF state is stabi- 
lized for 0.8 < n < 0.9, where AE is small. Hence the 
phase separation arises. By contrast, the phase separa- 
tion seen for a weak U in the ladder in Figs. [3] (a) (b) has 
a simple origin, i.e., both majorities and minorities are 
accumulated by the trapping potential and form doubly 
occupied states because the inter-particle interaction is 
weak. Finally, the spin-imbalance enhances the forma- 
tion of the finite hole-density Nagaoka ferromagnetism, 
i.e., the FF regions with 0.8 < n < 0.9 appear when 
p = 0.6 and 0.8 [Figs. [Hc)(d)], although they disappear 
in the spin-balanced case [Fig. [2][a)] due to a small gap 
AE around n ~ 0.8 (Fig. [p [2^ 

Phase diagram — Having clarified the mechanism for 
the phase-separated magnetism, we now explore a phase 
diagram against U, V, and i±. First, we focus on the 
dependence on V and U for p = 0.8 with tx = 1 in 
Figs.[5^a). We characterize the FF region by A'ff, which 
is the number of sites having rij = Sj within 10~ 2 normal- 
ized by the maximum value of this quantity. We find that 
there is optimum U ~ 50 and V ~ 3 to make N-p-p max- 
imum. To trace back how the optimum U and V arise, 
we can look at rii and Si plotted for various values of V 
and U in Figs. E^b)-(d). If we first combine Figs. [21(d), 
[5^b)(d) for the effect of varied V with a fixed U = 50, 
we can see that the central FF region broadens when V 
increases form 1 to 3 [Figs.[2td),[5Ub)] because the region 
with rii > 0.8 becomes wider with V. If the trap becomes 
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FIG. 5. (Color online) (a) Nff plotted against V and U for 
t j_ = 1 and p — 0.8. In the color-coding JVff is normalized by 
the maximum [a star at V = 3, U — 50, which corresponds to 
Fig.^d)] in the parameter region considered. Also plotted are 
Hi (blue solid lines) and s< (red dashed) with (V", U) = (1, 50) 
(b), (5,50) (c), and (3,10) (d). 
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FIG. 6. (Color online) (a) Color-coded Nff against V and tj 1 
for U = 50 and p = 0.8. Also plotted are m (blue solid lines) 
and St (red dashed) with (V",^ 1 ) = (1.5,0.6) (b), (1.5,0.2) 
(c), (2.5,0.2) (d). 



too strong in Fig.[5jd), however, the FF region gives way 
to the PM Mott plateau because ni reaches 1. If we turn 



to dependence on the interaction strength, we can look at 
the results for U = 10-50 with a fixed V = 3 [Figs. [2fd), 
(He)]. As the interaction is decreased to 10 in Fig. [5jc) 
the FF region shrinks, where the reason should be that 
AE becomes large at n ~ 0.8 when U is as small as 10. 

Finally we consider the effect of i±_ . A phase diagram 
in terms of Nff, against V and tj 1 this time, is displayed 
in Fig.[6fa). We find that a finite t±_ is suitable for largest 
A^ff in this figure, which is contrast to the homogeneous 
system in which the infinite i± is best. The reason is 
that a charge gap opens at = 0.5 [13], at which the FF 
phase gives way to the insulating state when ij 1 is small 
[Figs. HJfb)-(d)]. This degrades the FF phase stability. 

Summary — We have found for spin-imbalanced and 
strongly repulsive-interacting Fermi atoms on an opti- 
cal ladder that ferromagnetic regions appear in a phase- 
separated structures. Their emergence is caused by a 
combined effect of the strong interaction and the con- 
nectivity condition, i.e., the extended Nagaoka mecha- 
nism. The origin of the phase separation is explained by 
the density dependence of the energy gap between fully 
spin-polarized and other states in homogeneous system. 
As for finite temperatures, it was shown, with the dy- 
namical mean-field theory, that the finite hole-density 
Nagaoka ferromagnetism is robust against thermal fluc- 
tuations in two-dimensional system ;2G|. Then we ex- 
pect the phase separation is also robust at finite tem- 
peratures. Another future problem is how the itinerant 
ferromagnetism treated here in the spin-imbalanced con- 
dition would be related to the balanced case. Along the 
line, a study of correlation functions in spin-balanced case 
is under way 25j , where we observe phase separation be- 
tween non-magnetic and ferromagnetic (for the S x and 
S y direction) phases. We can however emphasize that 
the phase separation involving the FF regions for the S z 
direction will be easier to be observed with in situ imag- 
ing methods [l4| than the correlation functions. 
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